Multi-angle effects in self-induced oscillations for different supernova neutrino fluxes 
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The non-isotropic nature of the neutrino emission from a supernova (SN) core might potentially 
affect the flavor evolution of the neutrino ensemble, via neutrino-neutrino interactions in the deepest 
SN regions. We investigate the dependence of these "multi-angle effects" on the original SN neutrino 
fluxes in a three-flavor framework. We show that the pattern of the spectral crossings (energies 
where F Ve — F Vx , and F- e = F- x ) is crucial in determining the impact of multi-angle effects on 
the flavor evolution. For neutrino spectra presenting only a single-crossing, synchronization of 
different angular modes prevails over multi-angle effects, producing the known "quasi single-angle" 
evolution. Conversely, in the presence of spectra with multiple crossing energies, synchronization is 
not stable. In this situation, multi-angle effects would produce a sizable delay in the onset of the 
flavor conversions, as recently observed. We show that, due to the only partial adiabaticity of the 
evolution at large radii, the multi-angle suppression can be so strong to dramatically affect the final 
oscillated neutrino spectra. In particular three-flavor effects, associated with the solar parameters, 
could be washed-out in multi-angle simulations. 

PACS numbers: 14.60.Pq, 97.60.Bw 
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I. INTRODUCTION 

The characterization of the flavor conversions for 
neutrinos emitted by a stellar collapse is a field of 
intense activity. In particular, the flavor transformation 
probabilities in supernovae (SNe) not only depend on 
the matter background [l], , but also on the neutrino 
fluxes themselves: neutrino-neutrino interactions provide 
a nonlinear term in the equations of motio n KM i that 
causes collective flavor transformations ja-llq. Only 
recently fl7l - fli| it has been fully appreciated that in 
the SN context these collective effects give rise to 
qualitatively new phenomena (see, e.g., [20l ] for a recent 
review). The main consequence of this unusual type 
of flavor transitions is an exchange of the spectrum 
of the electron species v e {u e ) with the non-electron 
ones v x (p x ) in certain energy intervals. These flavor 
exchanges are called "swaps" marked by the "splits", 
which are the boundary features at the edges of each 
swap interval fill , l20| - [3l| . The location and the number 
of these splits, as well as their dependence on the 
neutrino mass hierarchy, is crucially dependent on the 
flux ordering among different neutrino species [HI, l33j . 

In this context, one of the main complication in the 
simulation of the flavor evolution is that the flux of 
neutrinos emitted from a supernova core is far from 
isotropic. The current-current nature of the weak- 
interaction Hamiltonian implies that the interaction 
energy between neutrinos of momenta p and q is 
proportional to (1 — v p ■ v q ), where v p is the neutrino 
velocity 0, |34| . In a non-isotropic medium this velocity- 
dependent term would not average to zero, producing 
a different refractive index for neutrinos propagating on 
different trajectories. This is the origin of the so-called 
"multi-angle effects" [lH , which hinder the maintenance 
of the coherent oscillation behavior for different neutrino 
modes [2l], I35rl37| . In [H| it has been shown that in 



a dense neutrino gas initially composed of only v e and 
V e with equal fluxes, multi-angles effects would rapidly 
lead to flavor decohcrence, resulting in flux equilibration 
among electron and non-electron (anti)neutrino species. 
On the other hand, in the presence of relevant flavor 
asymmetries between v e and v e multi-angle effects can 
be suppressed. In particular, during the early SN 
accretion phase, one expects as neutrino flux ordering 
®u e > $v e > = $ T7 X (38HIoj, defined in terms 
of the total neutrino number fluxes $^ for the different 
flavors. This case would practically correspond to 
neutrino spectra with a single crossing point at E — > oo 
(where F Ve = F Ux , and F„ e = F^ a ) since F Ue (E) > 
F Ux (E) for all the relevant energies (and analogously 
for v). For spectra with a single crossing, it has been 
shown in [36| that the presence of significant flavor 
asymmetries between v e and v e fluxes guarantees the 
synchronization of different angular modes at low-radii 
(r < 100 km), so that essentially nothing happens close 
to the neutrinospherc because the in-medium mixing 
is extremely small. Therefore, the possible onset of 
multi-angle effects is delayed after the synchronization 
phase. Then, the flavor evolution is adiabatic to produce 
spectral splits and swaps but not enough to allow the 
multi-angle instability to grow and produce significant 
decoherence effects. Therefore, the resultant neutrino 
flavor conversions would be described by an effective 
"quasi single-angle" behavior. In this case, the self- 
induced spectral swaps and splits would be only marginal 
smeared by multi-angle effects, as explicitly shown in (lH . 
l2ll [22J . This reassuring result has been taken as granted 
in most of the further studies that typically adopted the 
averaged single-angle approximation. 



However, this nice picture does not represent the 
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end of the story for multi-angle effects. 1 A different 
result has been recently shown in [42| |. where multi- 
angle effects are explored, assuming a flux ordering of the 
type § Vx > § Ve > $i7 e , possible during the SN cooling 
phase, where one expects a moderate flavor hierarchy 
among different species and a "cross-over" among non- 
electron and electron species is possible [381 - 4401 ] . This 
case would correspond to neutrino spectra with multiple 
crossing points, i.e. with F Ve (E) > F Ut (E) at lower 
energies, and F Vi! {E) < F V:i (E) at higher energies (and 
an alog ously for v). For such spectra, it has been shown 
in [43f that the synchronization is not a stable solution 
for a neutrino gas in presence of a large neutrino density. 
Therefore, collective flavor conversions would be possible 
at low-radii in the single-angle scheme 43[, in a region 
where one would have naively expected synchronization. 
However, it has been shown in [42| that the presence 
of a large dispersion in the neutrino-neutrino refractive 
index, induced by multi-angle effects, seems to block the 
development of these collective flavor conversions close 
to the neutrinosphcrc. This recent result extends the 
finding obtained with a toy model in [43j. The delay of 
the self-induced flavor conversions for this case is also 
visible in the multi-angle simulations in [44| . So, it is 
apparent that multi-angle effects are relevant not only for 
fluxes with small flavor asymmetries, where they trigger 
a quick flavor decoherence, but also in cases of spectra 
with multiple crossing points, where multi-angle effects 
can suppress flavor conversions at low-radii. 

Triggered by the contrasting impact of the multi-angle 
effects for fluxes typical of the accretion and cooling 
phase, we take a closer look at the dependence of these 
effects on the neutrino flux ordering. The plan of 
our work is as follows. In Sec. II we introduce our 
supernova flux models, and describe the equations for 
the flavor conversions in the multi-angle and single- 
angle case. In Sec. Ill we show and explain our 
numerical results for the single-angle and multi-angle 
flavor evolution for some representative choices of SN 
neutrino fluxes. In particular, we select three cases 
corresponding respectively to a) single-crossed neutrino 
spectrum with 3> 3>j7 e S> i producing a 

"quasi single-angle" flavor evolution, b) multiple-crossed 
spectrum with $„ T > $„ e > $y7 c , where single-angle 
and multi-angle evolutions give significantly different 
final neutrino spectra, c) small flavor asymmetries, i.e. 
5?t/ e ~ ~ , where the multi-angle suppression 
is small, and multi-angle decoherence produces a partial 
flavor equilibration among the different species. Finally, 
in Sec. IV we draw inferences from our results and 
summarize. Technical aspects are discussed in the 
Appendix. 



1 We mention that recently multi-angle effects have been included 
also in the study of the flavor evolution of the u e neutronization 
burst in O-Ne-Mg supernovae. We address the interested reader 

to El. 



II. SUPERNOVA NEUTRINO FLUXES AND 
EQUATIONS OF MOTION 

A. Supernova flux models 

In the presence of neutrino-neutrino interactions the 
flavor conversions for SN neutrinos are described by non- 
linear equations. Therefore, SN neutrino oscillations 
will crucially depend on the initial neutrino fluxes. We 
define F Va as the number flux of a given neutrino 
species v a emitted with energy E in any direction at 
the neutrinosphere. In a supernova v e and v e are 
distinguished from other flavors due to their charged- 
current interactions. The v lll v r and their antiparticles, 
on the other hand, are produced at practically identical 
rates. Following the standard terminology [45| . we define 
the two relevant non-electron flavor states as v XtV = 
cos^23^Tcos023' / r, where #23 — 7r/4 is the atmospheric 
mixing angle. Since the initial v x and v y fluxes are 
identical, the primary neutrino fluxes are best expressed 
in terms of v ei v e and v x . For half-isotropic emission 
these three relevant SN v original neutrino number fluxes 
for the different species are given by 



where 



4tt 2 R 2 {E v 



(1) 



(2) 



is the total number flux at the neutrinosphere radius R, 
that in our numerical examples we will take R = 10 km. 
The neutrino luminosity is L Ua and the neutrino average 
energy (E Va ) . The function tp Va (E) is the normalized 
neutrino spectrum (J dE ip Va — 1) and parametrized 
as 0: 



-PE/(E„ a ) 



(3) 



where j3 is a spectral parameter, and r(/3) is the Eulcr 
gamma function. The values of the parameters are model 
dependent (e.g. see the Fig. 3 in [32(). For our numerical 
illustrations, we choose 

«£„.)> ( e pJ) = ( 12 > 15, 18) MeV , (4) 

and P = 4, from the admissible parameter ranges fiij . 

We will consider three representative cases for the 
ratios of the neutrino fluxes, namely 



(5) 



The first case represents a flux ordering with a u e 
dominance, typical of the accretion phase, practically 
producing a single-crossed spectrum. The other two cases 
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represent fluxes possible during the cooling phase, with a 
moderate v x dominance and different flavor asymmetries, 
producing a multiple-crossed spectrum. We will see how 
multi-angle effects would have a different impact for these 
three cases. 



B. Equations of motion 

Mixed neutrinos are described by matrices of density 
p p and p p for each (anti)ncutrino mode. The diagonal 
entries are the usual occupation numbers whereas the 
off-diagonal terms encode phase information. We are 
studying the spatial evolution of the neutrino fluxes in 
a quasi-stationary situation. Therefore, the matrices p p 
do not explicitly depend on time, so that the evolution 
reduces to the Liouville term involving only spatial 
derivatives. Moreover, we assume spherical symmetry 
so that the only relevant spatial variable is the radial 
coordinate r. In this case, the explicit form of the 
equations of motion (EoMs) has been obtained in [H, H, 

GS H, S3, m. 



iv p • V r p p = [Hp,p p ] , (6) 

where v p is the velocity and the Hamiltonian reads 

Hp = ft p + Vmsw + V„„ . (7) 

In a three flavor scenario, the matrix of the vac- 
uum oscillation frequencies for neutrinos is ft p = 
diag(TO 2 , to 2 ,, m|)/2|p| in the mass basis. For antineu- 
trinos ft p — > — ft p . It will prove convenient to cast the 
matrix of vacuum oscillation frequencies in its traceless 
form 



8 — a— A3 



y/3 



(8) 



where A3 and Ag arc the two diagonal Gell-Mann 
matrices, which read respectively (45j 

A 3 = diag(l,-l,0) , 

A 8 = -^diag(l, 1, -2) . 

The vacuum oscillation frequency 

(9) 



2E 



is associated to the atmospheric mass-square difference 



Am 



to J 2 - The mass hierarchy parameter is 



a = sgn(ATO 2 tm ) 



A? 



(10) 



Ato 2 o1 — m\ — to 2 being the solar mass-square difference. 
Am^ tm > defines the normal mass hierarchy (NH), 
while Am^tm < the inverted hierarchy (IH). For 



the numerical illustrations, we take the neutrino mass- 
squared differences to be |Am^ tl J = 2 x 10~ 3 eV 2 and 
A' 



■in 



sol 



8 x 10 -5 eV 2 , close to their current best- 
fit values [49}. In SN neutrino flavor conversions, the 
parameters (Am 2 tm ,#i3) are responsible of conversions 



between 



y states, while 



(Ato 2 o1 ,#i 2 ) determine 



conversions between e — x states [451 ] . The mass hierarchy 
implies the dominance of e — y conversion effects over the 
e — x ones. However, are possible situations in which one 
finds an interesting in terp lay between the "atmospheric" 
and "solar" sectors [30, |31| . as we will show in the 
following. 

The matter effect in Eq. ([7]), due to the background 
electron density n e , in the weak interaction basis 
(y e , v^,v T ) is represented by [l[ 



V 



MSW 



v / 2G , Fn e diag(l,0,0) 



(11) 



Except at very early times (t < 300 ms) when the 
effective electron density n e would become larger than 
the neutrino density n v suppressing the self-induced 
flavor conversions [50j , one can account for matter effects 
in the region of collective oscillations just by choosing 
small (matter suppressed) mixing angles [17], Ell' which 
we take to be 6*13 = 612 = 10~ 3 . Matter effects in the 
region of collective oscillations (up to a few 100 km) 
also slightly modify the neutrino mass-square differences. 
Therefore, we take the effective mass-square differences 



A? 



= Am ttm COS 6»i3 , 



Am 2 ol cos 6<i2 



0.4Ar 



sol 



'sol 



[26|, l50|. Mikheyev-Smirnov- 



Wolfenstein (MSW) conversions typically occur after 
collective effects have ceased [U, Hf| ■ Their effects then 
factorize and can be included separately @]. Therefore, 
we neglect them in the following. 

Finally, the effective potential due to the neutrino- 
neutrino interactions is given by H3| 



V„„ = V2G F J 7^31 (Pq - Pq)(l - V F 



Vq) , (12) 



where the factor (1 — v p ■ v q ) implies multi- angle effects 
for neutrinos moving on different trajectories [HJ], as 
explained in Sec. I. 

We will solve the EoM's in the multi-angle case and 
compare the behavior of the flavor evolution with the 
solution obtained in the single- angle approximation [l8| . 
This latter requires the occurrence of the self-maintained 
coherence of the neutrino ensemble, i.e. that at a given 
location all the neutrino modes are aligned with each 
other, assuming they were aligned at the source [47| . 
In this case one obtains as angle-averaged EoM for 
the different energy modes, classified in terms of the 
frequency w 



ldrPoj = [H u ,p^ 



(13) 



where we have defined the "reduced neutrino density 
matrix" as [2Q| 



-p q if uj > 
if to < , 



(14) 
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whose diagonal components 



F Va {uj,r) 



(15) 



are normalized to the sum of the fluxes of all the neutrino 
species F v (ui, r) = F Ve (E, r) + F Vx (E, r) + F Vy (oj, r) (and 
analogously for antineutrinos). We clarify that when we 
use the w-variable, with the notation F v (oj) we would 
mean F U (E) x dE/doo. 

Neglecting the matter effects, the single-angle Hamil- 
tonian reads 143 



(16) 



In this equation the vacuum oscillation frequency oj when 
projected over the radial direction becomes J4?| 



OJ 
(V r ) 



where 



M = - 2 




(17) 



(18) 



is the angle-averaged radial neutrino velocity. The 
modification of the vacuum oscillation frequencies is 
relevant only near the ncutrinosphcre, therefore in the 
following we neglect it in our estimations. 

The neutrino-neutrino interaction term depends on the 
matrix of total density 



1 



*S + J- 



dojF u p u 



(19) 



normalized to the sum of the total neutrinos and 
antineutrinos flux at the neutrinosphere. 

The radial dependence of the neutrino-neutrino inter- 
action strength can be written as [2l|, H3] 



R2 r 



where 



(20) 



(21) 



represents the strength of the neutrino-neutrino poten- 
tial at the neutrinosphere. The r~ 2 scaling comes from 
the geometrical flux dilution, and the collinearity factor 



Cr 



1-y/l- [R/rf 



[R/rf 



for r 



(22) 



arises from the (1 — cos (9) structure of the neutrino- 
neutrino interaction. The asymptotic behavior for large 



r agrees with what one obtains by considering that all 
neutrinos are launched at 45° to the radial-direction (3|| . 
The known decline of the neutrino-neutrino interaction 
strength, fj,* ~ r~ 4 for r >• R, is evident Q . 

The behavior of the neutrino ensemble at large 
densities in the single-angle case can be characterized 
in terms of the invariants of the system. Namely, two 
lepton-numbers, which for small in-medium mixing angles 
read in flavor basis [2a, u5l| 



C 3 = 
Cr = 



1 



$0 i <J>Q 

U 1 V J — oo 



dojF v Tr(p u X a ) 



1 f +aa 
—— -o / dwF.Tr^Ag) , 



and the effective energy of the system [35l 



where 



1 



dcu oj F v p u 



^ V 1 ^ v J —OO 

is the first momentum of the density matrix. 



(23) 



(24) 



(25) 



III. MULTI-ANGLE EFFECTS FOR 
DIFFERENT NEUTRINO FLUXES 

In this Section we compare the results of the SN 
neutrino flavor evolution, obtained using the single-angle 
[Eq. (13])] and multi-angle [Eq. ©] EoM's for the three 
representative SN flux models introduced in Eq. ([5])- For 
the sake of brevity, we will show only results of the 
supernova neutrino flavor evolution in the case of inverted 
neutrino mass hierarchy. We think that this case is more 
interesting, since possible three-flavor effects, associated 
with Am^i have been recently found in the single-angle 
scheme for SN neutrino fluxes with multiple crossing 
points [13, [UJ. We will show that multi-angle effects 
can strongly suppress the impact of these three-flavor 
conversions. 

In our numerical simulations we fix the value of 
the neutrino-neutrino interaction strength at the neutri- 
nosphere 



fi*(R) = 2.1 x 10 6 km" 



(26) 



unless otherwise stated. This choice will imply that 
neutrino luminosities in the three different cases will be 
(in units of 10 51 erg/s) 

L Ve = 2.40 , = 2.00 , L Vx = 1.50 , 
L Ve = 1.20 , L Pe = 1.34 , L Vx = 2.14 , 
L Ve = 1.16 , L 9e = 1.41 , L Vx - 2.14 . 

(27) 
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FIG. 1: Case with $° e : $? e : <E>° x = 2.40 : 1.60 : 1.0. 
Three-flavor evolution in inverted mass hierarchy for the 
single-angle case for neutrinos (left panels) and antineutrinos 
(right panels). Upper panels: Initial energy spectra for v e 
(long-dashed curve) and u x ,y (short-dashed curve) and for 
v e after collective oscillations (solid curve). Lower panels: 
probabilities P ee (solid red curve), P ey (dashed blue curve), 
P ey (dotted black curve). 



In order to get stable results in our multi-angle simula- 
tions we typically use ~ 1500 angular modes. Finally, to 
saturate self-induced oscillation effects, we integrate the 
EoM's till r = 10 3 km. Technical details are discussed 
in the Appendix: in the following we focus only on the 
results and their interpretation. 



A. Spectrum with a single crossing 

We start our investigation with the SN neutrino flux 
ordering $° e : : $^ = 2.40 : 1.60 : 1.0, as 

representative of the accretion phase. This case has 
been the benchmark for most of the previous multi- 
angle studies (see, e.g., [HI, HH). As known, in this 
case the dynamics can be studied into the e — y two- 
neutrino system associated with (Am^ tm ,0i3). Three- 
flavor effects in the e — x sector, associated with 
(Am s 2 ol ,6»i 2 ) are negligible [Mill. 

In Fig. 1 we represent the initial neutrino fluxes at the 
neutrinosphcrc for all the different species and the final 
electron (anti)neutrino fluxes after collective oscillations 
(at r = 10 3 km) in the single- angle approximation (upper 



FIG. 2: As in Fig. 1, but for the multi-angle case. 



panels) and the corresponding P ee = P(v e — > Ve), 
P ex = P{y e — > v x ) and P ey = P(y e — > v y ) conversion 
probabilities (lower panels). Neutrinos are shown in the 
left panels and antineutrinos in the right panels. The 
corresponging results for the multi-angle case are given in 
Fig. 2. We stress that, except in the high-energy tails, the 
original neutrino spectra in the case under study always 
have an excess of electron (anti)neutrinos over the non- 
electron species, i.e. F Vll > F„ x and Fp e > F Vi . In the 
frequency variable — oo < uj < +oo, the crossing point at 
E -> oo (a; = 0) and the other two at E > 20 MeV 
appear so narrowly spaced, that the neutrino spectra 
superficially appear with only a single crossing point at 
id = 0, where F Ue = F^ = and F Vc =Fu x = 0. We 
address the interested reader to Fig. 3 of [231 and to the 
relative discussion on the case neutrino spectrum effectly 
behaving as a single-crossing one. We will show how this 
property is crucial in characterizing the flavor evolution. 

Concerning the neutrino fluxes after self-induced 
oscillations, in the single-angle case one finds a swap 
between v e and v y spectra above E ~ 10 MeV, producing 
the typical split in the final neutrino spectra. For 
antineutrinos, the swap between v e and v y is almost 
complete, the splitting energy E ~ 2 MeV being very low. 
Neglecting this low-energy feature [22| , the position of the 
v e split can be calculated using the conservation of the 
lepton number £ 8 in Eq. ([23]) (see, e.g., [H-il]). In the 
multi-angle case, the swap features remain unchanged, 
except for the smearing of the low-energy anti-neutrino 
spectral split, shown also in (22j . 



6 



neutrinos 



antineutrinos 



neutrinos 



antineutrinos 






^ 1 


= 1 


0.91 MeV- 




. 11 MeV 

- 43 MeV 






5 MeV 

i 



10' 
r(km) 



10 J 






I 1 






L 0.91 MeV - 


: 




/C 11MeV " 


5 MeV " 

i 


_ 43 MeV 




10* 
r(km) 



10 J 



FIG. 3: Case with : $g c : $° x = 2.40 : 1.60 : 1.0. Three- 
flavor evolution in inverted mass hierarchy in the single-angle 
case. Radial evolution of the diagonal components of the 
density matrix p for neutrinos (left panels) and antineutrinos 
(right panels) for different energy modes. 



In Fig. 3 and 4 we represent the radial evolution 
of the diagonal elements of the density matrix p ee , 
Pyy, Pxx for different energy modes for neutrinos (left 
panels) and antineutrinos (right panels). In particular, 
in the p ee panels the order of the energy modes is 
E = 0.91,5,11,43 MeV going from the curve starting 
with the highest value to the lowest one. This order is 
reversed in the p yy and p xx panels. Figure 3 represents 
the single-angle evolution, while Fig. 4 is for the multi- 
angle evolution, where the density matrix elements have 
been integrated over the angular distribution. Except 
for very low-energy antineutrino modes, we find that the 
evolution of the density matrix is rather similar in the 
single-angle and multi-angle case. We find the presence 
of synchronized oscillations 0, [lj| till r ~ 85 km. Till 
there all the p u stay pinned to their original value. 

Significant flavor conversions start only after synchro- 
nization when bipolar oscillations fl9j start to swap the 
flavor content of the system. The synchronization radius 
can be found through the condition [2l[ 



Pr > 



4 uj s Tr(aA 8 ) 
V3 [Tr(p\sW 



(28) 



Using the alignment approximation for the two blocks 
of neutrinos and antineutrinos, the synchronization 



FIG. 4: As in Fig. 3, but for the multi-angle case. 

frequency for e — y conversions, reads pll ] 
J+°°<h> q;(F°(£)-F°(a,)) 

f+°°du, W (F°>)-F°>)) 



2Ar 



2(*g.-*g B ) 
L v J(E v .)*-L u J(E Vm ? 



L„J{E„ e )-L v J(E u J 



I^JiE^y-L^/jE^ 
LuJ(Ev e ) -L V J(E V J 
~ 0.68 km" 1 , 

for our input fluxes. The function 

1 



$0 i <J>0 



du sgn(u>)F u p u 



Then 



Tr(<xA 8 ) 



Pee H~ Pee 

V3 



1 L Ve (E Uc )(E Vx )+ Lj7 e )(E Vx )-2L Vx (E Uc )(E Ve ) 
v/3 L Ve {E Vx }(E Ve } + L Ve {E Vc ) {E Vx ) + 4L„ x {E Ve ) (£^ e ) 
0.144 , 



and 



Tr(pA 8 ) 



Pee Pee 

~7T~ 
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1 



(E v „)[L v SEv c ) - Lv c (E v .)\ 



V3 L Ve (E Vx }(E Ve } + {E Ve ) (E„ x ) + AL Ux (E Ve ) (E Ve 
0.06 . 



These numbers imply that bipolar conversions would 
start when p* ~ 100 uih, i-e. at r ~ 85 km, as observed 
in the Fig. 3. 

We stress that the presence of synchronization at low- 
radii in this case is crucially related with the original 
spectrum with a single crossing point. As discussed 
in [43l | . energy conservation at large neutrino densities 
(i.e. when p* — > 00) implies that p 2 and thus p 
are conserved [see Eq. ([24]) ]. and therefore behave as 
collective objects. If p is conserved, energy conservation 
implies that also Tr(pX s ) has to be conserved (neglecting 
the subleading three-flavor effects associated to Am^). 
For a single-crossed spectrum, the quantity Tv(pXs) — 

[(pee — Pyy) + {p ee ~ Pyy)]/V^ is maximal [29|]. Therefore 
moving any energy mode p^ from its initial value 
would make the energy of the system larger. In this 
situation, each p^ must remain pinned with the global 
density matrix p. Therefore, only the synchronization 
among different modes is possible [43||. Since during the 
synchronization phase the different p u remain aligned to 
their original value, multi-angle affects are suppressed. 

Then, during the phase of the spectral swapping, the 
spectrum near the crossing point acts like an inverted 
pendulum J2||. The swap sweeps through the spectrum 
on each side of the crossing, and the modes at the edge 
of the swap precess at an average oscillation frequency 
k ~ Arng tm /2E ~ 0.5 km -1 for a typical energy 
E ~ 10 MeV in the region of the swap. Since, the 
length scale 1^ = \d\x\ p*{r) / dr]' 1 = r/4 is larger than 
the evolution is adiabatic concerning the swapping 
dynamics [13, [H|. However, the time-scale of multi-angle 
effects is determined by the bipolar oscillation frequency 
wh ~ y/2ufi* ~ r~ 2 [35} that decreases faster than 
l^ 1 . Qualitatively, the effect of the synchronization is 
to postpone conversions at large radii, where multi-angle 
effects would require more "time" to develop since p* 
is smaller, and the relatively fast decrease of p* would 
reduce the adiabaticity for multi-angle effects. As a 
consequence these do not grow significantly This is 
qualitatively the origin of the so called "quasi single- 
angle" behavior (36| . observed for neutrino flux ordering 
typical of the accretion phase. We have explicitly 
checked that modifying artificially the adiabaticity of the 
evolution, i.e. choosing a very slow decrease of p* multi- 
angle decoherence is unavoidable also in the case we are 
considering. Therefore, as suggested in (36|, the absence 
of multi-angle decoherence during the accretion phase 
seems due to a luckily conspiracy of different time-scales. 
Unfortunately, till now it has not been developed yet a 
complete theory for that. Therefore, the interpretation 
of these effects is mostly based on the experience gained 
through numerical observations. 



B. Spectrum with multiple crossings 

We pass now to the case <&° c : $" e : $° x = 0.85 : 
0.75 : 1.0 representative of recent simulation results for 
the neutrino flux ordering during the cooling phase [38| . 
This case has been widely studied in the single-angle 
approximation in j29U3"i| and corresponds to spectra that 
in uj variable present well separated multiple crossings, 
in which multiple spectral splits can arise around the 
crossing points (see Fig. 2 of [29[). Moreover, in this case 
peculiar three-flavor effects, associated with conversions 
between e and x states have been recently discussed (30l . 
[3l| . Multi-angle simulations in this case have shown 
a smearing of the splitting features [29|, E2I and a 
suppression the flavor evolution at low-radii |42|, EJ • 

In Fig. 5 we show the initial (anti)neutrino spectra 
for all the flavors and the final electron (anti) neutrino 
ones for the single-angle case (upper panels) and the 
corresponding conversion probablities at the end of the 
flavor evolution. Multi-angle results are shown in Fig. 6. 
Starting with the single-angle case, we see that both 
v e <-> v y as well as v e f-> v y swaps appear at intermediate 
energies 5 MeV < E < 20 MeV. Moreover, for E > 25 
MeV, there are additional v e <-> v x and v e <-> v x swaps. 
As result of this complex dynamics, in the single-angle 
case the oscillated v e spectrum shows a single split, the 
one at low-energies, producing the swap with v y after 
the split, while the high-energy split is canceled by the 
further swap between v & and v x . We observe that in the 
antincutrino sector the spectral swaps are not complete, 
due to smaller spectral differences between v e and v x that 
lead to a less adiabatic evolution, as explained in [3l| . 

Passing to the multi-angle case, the differences with 
the single-angle case in the final neutrino spectra and in 
the conversion probabilities are significant. In general, 
all the splitting features are less pronounced than in the 
single-angle case. Starting with the neutrino case, we 
see that the low-energy spectral split is smeared. More 
remarkably, the high-energy swap between v e and v x is 
strongly suppressed. As a consequence, the final neutrino 
spectra keep track only of the swap between v e and v yi 
which now is less sharp than in the single-angle case. This 
effect would produce a high-energy splitting feature in 
the final v e spectrum at E ~ 25 MeV, not present in the 
single-angle three-flavor calculations. For antineutrinos, 
the impact of multi-angle effects is even stronger: flavor 
transformations are significantly suppressed, except in 
the energy range close to crossing point of the spectra, 
where two-flavor D e — > D y transformations take place. 

To explain the above observations, we represent the 
radial evolution of the diagonal elements of the density 
matrix p ee , p yy , p xx for different energy modes, for 
neutrinos (left panels) and antineutrinos (right panels) in 
Figs. 7 and 8. Figure 7 represents the single- angle case, 
while Fig. 8 represents the multi-angle case. In this latter 
case, the density matrix elements have been integrated 
over the angular distribution. In the p ee panels the order 
of the energy modes is E = 3.3,5.7,19,30 MeV going 
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FIG. 5: Case with $° e : $g e : $° x = 0.85 : 0.75 : 1.0. 
Three-flavor evolution in the single-angle case for neutrinos 
(left panels) and antineutrinos (right panels). Upper panels: 
Initial energy spectra for v e (long-dashed curve) and u x , v 
(short-dashed curve) and for v e after collective oscillations 
(solid curve). Lower panels: probabilities P ee (solid red 
curve), P ey (dashed blue curve), P ey (dotted black curve). 



from the curve starting with the highest value to the 
lowest one. This order is reversed in the p yy and p xx 
panels. Starting with the single-angle case, we see that 
differently from the case studied in the Sec. Ill A, flavor 
conversions are possible at low radii (r > 30 km) in a 
region where we would have expected synchronization. 
The effect of the small in-medium mixing is only to 
logarithmically delay the onset of the flavor conversions. 

As explained in [43| . the crucial difference with 
respect to the previous case is that, applying the energy 
conservation [Eq. (|24p] to the case of multiple-crossed 
spectra, it is possible to conserve Tr(pAs) and aTr(/5A3), 
by flipping parts of the original neutrino spectra around 
the crossing points. Therefore, the synchronization 
behavior found in the case of single-crossed spectra is 



not necessarily stable in this case. Indeed, in 43 1 
it has been discussed the possibility of a novel form 
of flavor conversions for large p* in terms of a self- 
induced parametric resonance that would destabilize the 
synchronization. The origin of this effect is that since 
each p u would precess around the Hamiltonian H w with 
a frequency p*, itself must vibrate with the same 
frequency. In this situation, for large p* the total density 
matrix p continues to behave as a collective object, but 



FIG. 6: As in Fig. 5, but for the multi-angle case. 



is not static and vibrates itself with a frequency p*. For 
single-crossed spectra, the energy conservation prevents 
any seizable change in the different p w , therefore the 
synchronized behavior prevails. Conversely, in presence 
of multiple-crossed neutrino spectra, the different p^ 's do 
not remain aligned with the global p and start to librate 
relative to each other. 

We also note that in this case three-flavor effects, 
associated with (Am^ ol ,9i2) are present. In particular, 
from the figure, one realizes that collective librations first 
start in the e — y system (at r > 30 km), and then trigger 
the e — x conversions (at r > 60 km). The explanation of 
this dynamics has been recently given in [3l|, where we 
redirect the interested reader. 

Passing to the multi-angle case of Fig. 8, we find that 
now the low-radii flavor conversions are suppressed. This 
effect has been recently described in , even if there it 
has not made the connection between this behavior and 
the multi-angle suppression of the parametric resonance. 
This effect has been pointed out in [43[ for generic 
neutrino gases with half-isotropic neutrino distributions, 
in cases where in the single-angle scheme the parametric 
resonance was found. The point is that in the multi- 
angle case, for large p* the different p w 's would precess 
with different velocities whose spread, due to multi-angle 
effects, is much larger than w. As a consequence of 
this dispersion, there cannot exist a collective parametric 
resonance for all the neutrino modes. Therefore, the 
collective behavior found in the idealized single-angle 
case is "fragile" and easily suppressed by multi-angle 
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flavor evolution in inverted mass hierarchy in single-angle 
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(right panels) for different energy modes. 



effects in the realistic anisotropic supernova environment. 
In this long as multi-angle effects are dominant, 

they would suppress any flavor conversion at small radii. 
As discussed in [42j , the role of multi-angle effects in this 
case is analogous to the one that they have in presence 
of a large matter term (n e 3> n v ) |50[. There, multi- 
angle effects introduce a significant dispersion for the 
matter potential encountered by neutrinos on different 
trajectories, that once more would suppress self-induced 
conversions. 

In the presence of multi-angle suppression, we expect 
that flavor conversions would start only when neutrino- 
neutrino interactions are not strong enough to maintain 
the collective behavior of p [43 |. Neutrino modes would 
remain frozen to their initial condition, till the neutrino- 
neutrino interaction p* becomes comparable with the 
averaged vacuum oscillation frequency of the neutrino 
ensemble 1451 



FIG. 8: As in Fig. 7, but for the multi-angle case. 



Tr pA 8 ) 
Tr{p\ s ) 

where for our input spectra 
f>< 



Pt 



IPyy 



(29) 



Pee 2f)yy 



1.1 x IQ- Z km 



and 



Pee Pee 

L Ve (Ey. )(£„,) + Lg. (E Vc ) {E Vx ) - 2L„ x (Er U<i ) (E Vc ) 
L Ue (E Vi ) (E Ve ) + L Vc {E Vb ) (E Vm ) + AL Vm {E Vs ) (E Ve ) 
~ 1.6 x 10~ 2 , (30) 

which numerically leads to (cj) ~ 0.68 km -1 . 

We expect flavor conversions to start in the e — y sector 
when 



V3p*Tr(p\ s ) 



(31) 



LuJ{E Vc ) 



that would correspond to p* ~ 60 (lu). For our 
input spectra, we would get as onset radius r ons ~ 
92 km, in qualitative agreement with the simulation 
of Fig. 8. We note that numerically the conversions' 
onset radius in this case is very similar to the one of 
Sec. Ill A, even if the conditions that determine it are 
different, as well as the further flavor evolution. In 
particular, we expect that at large radii the evolution 
in this case would be less adiabatic than in the case of 
Sec. Ill A. As explained in [29|], in the final phases of 
the swapping dynamics, the neutrino and antineutrino 
spectra evolve quite independently, and the precession 
frequencies of the two blocks are not governed by a 
common p*, but by individual /i*'s proportional to 
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FIG. 9: Case with $° e : : = 0.81 : 0.79 : 1.0. 

Three-flavor evolution in the single-angle case for neutrinos 
(left panels) and antineutrinos (right panels). Upper panels: 
Initial energy spectra for v e (long-dashed curve) and u x , v 
(short-dashed curve) and for v e after collective oscillations 
(solid curve). Lower panels: probabilities P ee (solid red 
curve), P ey (dashed blue curve), P ey (dotted black curve). 



the flux differences \F„ e - F„J QF^ - F 9 J). They 
behave essentially as two uncoupled oscillators because 
the neutrino-neutrino interaction /i* is now smaller than 
the frequency difference of the two blocks. Since in this 
case the flux differences are smaller than in the previous 
example of Sec. Ill A, the effective precession frequencies 
are smaller, with a resultant less adiabatic evolution. 
This would explain the more pronounced smearing of the 
sharp splitting features observed in this case in the multi- 
angle simulations. 

The effect of the multi-angle delay is that flavor 
conversions start in a region in which the neutrino- 
neutrino interaction strength is weaker than in the 
corresponding single-angle case. The typical time-scale 
at which the off-diagonal components in the density 
matrix grow is given by the bipolar period which, for 
e — y transitions, is th ~ (2w^*) -1 / 2 The delay 

of the conversions means that they start at a smaller 
fi* where they require more time to develop. Moreover, 
going at larger r makes the evolution less adiabatic, 
resulting in weaker flavor conversions than in the single- 
angle case. This would explain the suppression of the 
e — x conversions. At this regard, we remind that 
e — x conversions are triggered by the dominant e — 



FIG. 10: As in Fig. 9, but for the multi-angle case. 



y conversions, via the #13 coupling between the two 
sectors |3lj . Due to the suppression of e — y conversions, 
the e — x transitions would also be delayed. The typical 
time-period at which the off-diagonal components in the 
e — x sector grow, is T£ ~ (aui^i*)^ 1 ^ 2 . The mass 
hierarchy implies that T£ ~ 8th- As a consequence of 
this slower growth of the e — x conversions, the multi- 
angle suppression is enough to make tl too slow to be 
effective. 

The suppression of conversions between v e and v x 
explains the appearance of the high-energy spectral 
split in v e observed in the multi-angle case of Fig. 7. 
Concerning antineutrinos, we realize that not only e — x, 
but also e — y conversions are strongly inhibited. We 
associate this behavior with the stronger violation of 
adiabaticity in the antineutrino sector, due to the smaller 
spectral differences among P e and v x [3lj . 

We stress that the adiabaticity plays a crucial role in 
determining the impact of the multi-angle suppression in 
the flavor evolution. Indeed, we explicitly checked that 
for the same flux ordering, increasing the adiabaticity, 
by increasing the neutrino-neutrino potential by a factor 
of five, the multi-angle suppression is less dramatic. In 
particular, e — x conversions reappear. Moreover, also 
the suppression of flavor oscillations in the antineutrino 
sector is relieved. We find that this result is in agreement 
with the case shown in [42j ■ There the choice of higher 
neutrino luminosities with respect to the ones we are 
considering as benchmark case, allows the multi-angle 
suppression to be strong only at low-radii, while the final 
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FIG. 11: Case with : $£ e : = 0.81 : 0.79 : 1.0. Three- 
flavor evolution in inverted mass hierarchy in single-angle 
case. Radial evolution of the diagonal components of the 
density matrix p for neutrinos (left panels) and antineutrinos 
(right panels) for different energy modes. 



spectra arc similar to the ones obtained in the single- 
angle case. 



C. Spectrum with small flavor asymmetries 

Finally, we consider the case : $? e : <t>° x = 0.81 : 
0.79 : 1.0 which is intended to represent a case which a 
flux ordering possible at late times, where asymmetries 
among v e and v e can become small [HI, H3| ■ The small 
asymmetry case has been pointed as representative of 
flavor decoherence associated to multi-angle effects [HI, 

m. 

In Fig. 9 show the initial (anti)neutrino spectra for 
all the flavors and the final electron (anti) neutrino 
ones for the single-angle case (upper panels) and the 
corresponding conversion probablitics at the end of 
the flavor evolution. Multi-angle results are shown 
in Fig. 10. The splitting features in the single-angle 
case are similar to the ones observed in Fig. 5 of Sec. 
Ill B. However, the multi-angle results are dramatically 
different. We observe that neutrino and antineutrino 
spectra tend toward flavor equilibration with different 
conversion probabilites displaced around 1/3. 

In Fig. 11 and 12 wc represent the radial evolution 
of the diagonal elements of the density matrix p ee , p yy , 



FIG. 12: As in Fig. 11, but for the multi-angle case. 



p xx for different energy modes for neutrinos (left panels) 
and antineutrinos (right panels). In the p ee panels the 
order of the energy modes is E = 2.5, 5, 19, 30 MeV 
going from the curve starting with the highest value 
to the lowest one. This order is reversed in the p yy 
and p xx panels. Figure 11 represents the single-angle 
case, while Fig. 12 represents the multi-angle case. The 
difference of the flavor evolution in the two cases is 
striking. In particular, multi-angle suppression blocks 
e-y flavor conversions till r ~ 40 km, as predictable from 
Eq. (|31|) . Then, since conversions start at small r where 
the evolution is more adiabatic than in the case of Sec. 
Ill B, also e-x oscillations have chance to develop at 
r > 60 km. The stronger adiabaticity in this case allow 
also multi-angle effects to have enough time to develop 
before the neutrino-neutrino interaction term becomes 
small. Then, the multi-angle effects smear the flavor 
conversions, producing a tendency towards a three-flavor 
decoherence of the ensemble in both v and 77 sectors. 



IV. CONCLUSIONS 

We have performed an exploration on the dependence 
of multi-angle effects in self-induced supernova neutrino 
oscillations on the ori gina l neutrino fluxes. Most of the 
previous studies [H, T2TI I22T I36T ) focused on neutrino 
fluxes typical of the accretion phase, with a pronounced 
v e excess, de facto behaving like spectra with a single 
crossing. In this situation, the synchronization of 
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TABLE I: Summary of multi-angle effects, 3v effects and spectral spits for different SN neutrino fluxes. 



Initial spectral pattern 


Multi-angle effects 


Amg 0l -effects 


Spectral splits 


Single crossing 


marginal 


absent 


robust 


Multiple crossings 


relevant 


present / absent 


smeared 


Small flavor asymmetries 


strong 


present 


washed-out 



different angular modes at low-radii prevails over multi- 
angle effects. Then, when flavor conversions start, these 
are adiabatic to produce the spectral swaps and splits, 
but not enough to allow multi-angle decoherence to 
emerge. The result is the known "quasi single-angle" 
evolution. 

However, one has to be cautious in generalizing this 
reassuring result. In this context, we have shown that 
multi-angle effects can produce significant deviations 
in the flavor evolution with respect to the three-flavor 
single-angle case, for neutrino fluxes with a moderate 
flavor hierarchy and a v x excess, as possbile during 
the supernova cooling phase. In this situation, the 
presence of multiple crossing points in the original 
neutrino spectra destabilizes the synchronization at large 
neutrino densities [43|. In absence of synchronization, 
in the single-angle scheme collective conversions would 
be possible at low-radii. However, multi-angle effects 
introduce a large dispersion in the neutrino-neutrino 
potential, that prevents any possible collective flavor 
conversion at low radii. As a consequence, in the 
multi-angle scheme there would be a significant delay 
of the onset of the flavor conversions, as recently 
observed [42j . We have shown that this multi-angle delay 
can produce dramatic changes, not only in the deepest 
supernova regions, as shown in |42| . but also in the final 
oscillated neutrino spectra. In particular, the multi- 
angle suppression can be so strong to allow the onset 
of the flavor evolution only at a large radius, when the 
evolution is less adiabatic. Depending on the violation of 
adiabaticity, in the inverted mass hierarchy there could 
be a suppression of the three-flavor effects, associated 
to the solar sector. This would dramatically change the 
pattern of swaps and splits in the final neutrino spectra. 
Finally, if the flavor asymmetries between u e and V e are 
very small, the multi-angle suppression occurs only close 
to the neutrinospherc. In this situation, the stronger 
adiabaticity of the evolution allows multi-angle effects 
to act efficiently also after the onset of the conversions, 
tending to establish a three-flavor equilibration in both 
neutrino and antincutrino sectors. 

In Table I we summarize our results on the role of 
multi-angle effects, 3^ effects and spectral splits for 
different SN neutrino fluxes. In this work we have 
explicitly shown numerical results only for the case 
of neutrino inverted mass hierarchy. However, we 
have checked that the impact of multi-angle effects is 
qualitatively similar also for the normal hierarchy case. 



From our numerical explorations, it results that self- 
induced flavor transformations of supernova neutrinos 
during the cooling phase arc a continuous source of 
surprises. The richness of the phenomenology, in the 
presence of neutrino spectra with multiple crossing 
points, was first realized in [2^] with the discovery of 
the possibility of multiple spectral splits in both the 
mass hierarchies. Then, it was realized that for these 
neutrino fluxes, three-flavor effects can play a significant 
role in inverted mass hierarchy, changing the splitting 
pattern expected from two- flavor calculations pOl . [3l| . 
Now, we show that also multi-angle effects are crucial 
in characterizing the flavor evolution in this case, 
and could potentially kill the three-flavor effects. In 
general, self-induced flavor conversions for spectra with 
multiple crossing points challenge most of the naive 
expectations on which was based the original picture of 
the collective supernova neutrino conversions: low-radii 
synchronization, subleading role of multi-angle and three- 
flavor effects. 

The discovery of these new effects adds additional 
layers of complications in the simulation of the flavor 
evolution for supernova neutrinos. In particular, our 
result shows that during the cooling phase three-flavor 
multi-angle simulations are crucial to obtain a correct 
result. Multi-angle effects would be taken into account 
to assess the impact of collective neutrino oscillations on 
the r-process nucleosynthesis in supernovae, as recently 
investigated in [53| . The impact of the multi-angle effects 
would crucially depend on different SN input: neutrino 
luminosities and flavor asymmetries, neutrinospherc 
radius, etc. Since all these quantities significantly change 
during the neutrino emission, one would expect time- 
dependent effects. At this regard, the possibility to 
detect signatures of these effects in the next galactic 
supernova neutrino burst (33| would motivate further 
analytical and numerical investigations. 



Appendix 

We discuss here a few technical aspects of the multi- 
angle numerical simulations, we performed on our local 
computer facility (with Fortran 77 codes running a 
Linux cluster with 48 processors per CPU with 128 
Gb of shared RAM memory). Equation ©, after 
discretization, provides a set of 16 x Ne x N u ordinary 
differential equations in r, where Ne and N u are the 
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number of points sampling the (anti) neutrino energy E 
and emission angle. In particular, we find convenient 
to label the neutrino angular modes in terms of the 
variable [36[ 

u = sm 2 e R , (32) 

where 8r is the zenith angle at the neutrino sphere r = R 
of a given mode relative to the radial direction. With 
this choice, the parameter u is fixed for every neutrino 
trajectory. 

We have then performed our simulations of the three- 
flavor neutrino evolution using a Runge-Kutta integra- 
tion routine taken from the CERNLIB libraries [54| . We 
fixed the numerical tolerance of the integrator at the level 
of 1CP 6 and increased the number of sampling points in 
angle and energy till we reach a stable numerical result. 
In this situation we estimate a numerical (fractional) 
accuracy of our results better than 10 ~ 2 . In order to 
have a clear energy resolution of the spectral splits we 
took Ne = 10 2 energy points, equally distributed in the 
range E £ [0.1, 80] MeV. The number of angular modes is 



also a crucial choice, since it is well known that a sparse 
sampling in angle can lead to numerical artifacts that 
would destroy the collective behavior of the neutrino self- 
induced conversions (36j . Typically, numerical stability 
would require N u = O(10 3 ) angular modes. 

Since we can claim to have reached stable numerical 
simulations, we are confident in the accuracy of the 
results obtained in this work. Moreover, we have been 
able to reproduce previous results presented in literature 
for cases similar to the ones we are investigating (see, 
e.g., PHI). 



Acknowledgments 

We thank Sovan Chakraborty, Basudeb Dasgupta 
and Georg Raffelt for reading the manuscript and for 
comments on it. The work of A.M. was supported 
by the German Science Foundation (DFG) within the 
Collaborative Research Center 676 "Particles, Strings 
and the Early Universe" . 



L. Wolfenstein, "Neutrino Oscillations In Matter," Phys. 
Rev. D 17, 2369 (1978); S. P. Mikheev and A. Yu. 
Smirnov, "Resonance Enhancement Of Oscillations In 
Matter And Solar Neutrino Spectroscopy," Yad. Fiz. 42, 
1441 (1985) [Sov. J. Nucl. Phys. 42, 913 (1985)]. 

A. S. Dighe and A. Y. Smirnov, "Identifying the neutrino 
mass spectrum from the neutrino burst from a super- 
nova," Phys. Rev. D 62, 033007 (2000) [hep-ph/9907423] . 
J. Pantaleone, "Neutrino oscillations at high densities," 
Phys. Lett. B 287, 128 (1992). 

G. Sigl and G. Raffelt, "General kinetic description of 
relativistic mixed neutrinos," Nucl. Phys. B 406, 423 
(1993). 

B. H. J. McKellar and M. J. Thomson, "Oscillating 
doublet neutrinos in the early universe," Phys. Rev. D 
49, 2710 (1994). 

Y. Z. Qian and G. M. Fuller, "Neutrino-neutrino 
scattering and matter enhanced neutrino flavor transfor- 
mation in Superno vae," Phys. Rev. D 51, 1479 (1995) 



astro-ph/9406073l 



Samuel, "Neutrino oscillations in dense neutrino gas- 
es," Phys. Rev. D 48, 1462 (1993). 

V. A. Kostelecky and S. Samuel, "Neutrino oscillations 
in the early universe with an inverted neutrino mass 
hierarchy," Phys. Lett. B 318, 127 (1993). 
V. A. Kostelecky and S. Samuel, "Self-maintained 
coherent oscillations in dense neutrino gases," Phys. Rev. 
D 52, 621 (1995) [hep-ph/9506262] . 

S. Samuel, "Bimodal coherence in dense selfinteract- 
ing neutrino gas es," Phys. Rev. D 53, 5382 (1996) 
[hep-ph/9604341] . 

ST Pastor, <G G. Raffelt and D. V. Semikoz, 
"Physics of synchronized neutrino oscillations caused 
by self-interactio ns," Phys. Rev. D 65, 053011 (2002) 
[hep-ph/0109035] . 

Y. Y. Y. Wong, "Analytical treatment of neutrino 



[13] 

[14] 

[15] 
[16] 

[17] 

[18] 

[19] 

[20] 
[21] 

[22] 



asymmetry equilibration from flavour oscillations in 
the early universe ," Phys. Rev. D 66, 025015 (2002) 
[hep-ph/ 0203180] , 

K. N. Abazajian, J. F. Beacom and N. F. Bell, "Stringent 
constraints on cosmological neutrino antineutrino asym- 
metries from synchronized flavor transformation," Phys. 



Rev. D 66, 013008 (2002) [arXiv:astro- ph/0203442|. 
S. Pastor and G. Raffelt, "Flavor oscillations in the 
supernova hot bubble region: Nonlinear effects of 
neutrino background," Phys. Rev. Lett. 89, 191101 
(2002) [astro-ph/0207281| . 

R. F. Sawyer, "Classical instabilities and quantum speed- 



up in the evolution of neutrino clouds," hep-ph/0408265 
R. F. Sawyer, "Speed-up of neutrino transformations in a 
supernova environ ment," Phys. Rev. D 72, 045003 (2005) 
[hep-ph/0503013] . 

H. Duan, G. M. Fuller and Y. Z. Qian, "Collective 
Neutrino Flavor Transformation In Supernovae," Phys. 
Rev. D 74, 123004 (2006) [astro-ph/0511275| . 
H. Duan, G. M. Fuller, J. Carlson and Y. Z. Qian, 
"Simulation of coherent non-linear neutrino flavor trans- 
formation in the supernova environment. I: Correlated 
neutrino trajectori es," Phys. Rev. D 74, 105014 (2006) 
[astro-ph/0606616l . 

S. Hannestad, G. G. Raffelt, G. Sigl and Y. Y. Y. Wong, 
"Self-induced conversion in dense neutrino gases: Pendu- 
lum in flavour space," Phys. Rev. D 74, 105010 (2006 ) 
[Erratum-ibid. D 76, 029901 (2007)] [astro-ph/0608695| . 

H. Duan, G. M. Full er and Y. Z. Q ian, "Collective 
Neutrino Oscillations." larXiv:1001. 27991 [hep-ph], 

G. L. Fogli, E. Lisi, A. Marrone and A. Mirizzi, 
"Collective neutrino flavor transitions in supernovae and 
the role of trajec tory averaging," JCAP 0712, 010 (2007) 
[arXiv:0707. 19981 [hep-ph]]. 

G. L. Fogli, E. Lisi, A. Marrone, A. Mirizzi and 

I. Tamborra, "Low-energy spectral features of supernova 



14 



(anti)neutrinos in inverted hierarchy," Phys. Rev. D 78, 

097301 (2008) [arXiv:0808.0"807l [hep-ph]] . 
[23] G. G. Raffelt and A. Y. Smirnov, "Self-induced spectral 

splits in supernova neutrino fluxes," Phys. Rev. D 76, 

081301 (2007) [ Erratum-ibid. D 77, 029903 (2008)] 

|arXiv:0705.1830l [hep-ph]] . 
[24] G. G. Raffelt, A. Y. Smirnov, "Adiabaticity and spectral 

splits in collective neutrino transformations," Phys. Rev. 

D76, 125008 (2007). | arXiv:0709.464T1 [hep-ph]] . 
[25] H. Duan, G. M. Fuller, J. Carlson et al., "Neutrino Mass 

Hierarchy and Stepwise Spectral Swapping of Supernova 

Neutrino Flavors ," Phys. Rev. Lett. 99, 241802 (2007). 

[arXiv: 0707.02901 [astro-ph]] . 
[26] H. Duan, G. M. Fuller, Y. -Z. Qian, "Stepwise spectral 

swapping with three neutrino flavors," Phys. Rev. D77, 

085016 (2008). larXiv:0801.f363l [hep-ph]]. 
[27] J. Gava and C. Volpe, "Collective neutrinos oscillation 

in mat ter and CP-violat ion," Phys. Rev. D 78, 083007 

(2008) |arXiv:0807.34T8l [astro-ph]] . 

[28] J. Gava, J. Kneller, C. Volpe and G. C. McLaugh- 
lin, "A dynamical collective calculation of supernova 
neutrino signals, " Phys. Rev. Lett. 103, 071101 (2009) 
larXiv:0902.03l7l [hep-ph]] . 

[29] B. Dasgupta, A. Dighe, G. G. Raffelt and A. Y. Smirnov, 
"Multiple Spectral Splits of Supernova Neutrinos," Phys. 
Rev. Lett. 103, 051105 (2009) |arXiv:0904.3542l [hep- 
ph]]. 

[30] A. Friedland, "Self-refraction of supernova neutrinos: 
mixed spectra and three-flavor instabilities," Phys. Rev. 
Lett. 104, 191102 (2010) |arXiv:1001. 09961 [hep-ph]]. 

[31] B. Dasgupta, A. Mirizzi, 1. Tamborra and R. Tomas, 
"Neutrino mass hierarchy and three-flavor spectral splits 
of supernova neut rinos," Phys. Rev. D 81, 093008 (2010) 
|arXiv:1002.2"943l [hep-ph]] . 

[32] G. Fogli, E. Lisi, A. Marrone and 1. Tamborra, "Su- 
pernova neutrinos and antineutrinos: ternary luminosity 
diagram and spectral split patterns," JCAP 0910, 002 

(2009) |arXiv:0907.5TT5l [hep-ph]] . 

[33] S. Choubey, B. Dasgupta, A. Dighe and A. Mirizzi, 
"Signatures of collective and matter effects on supernova 
neutrinos at large detectors," arXiv: 1008.0308 [hep-ph]. 

[34] J. T. Pantaleone, "Dirac neutrinos in dense matter," 
Phys. Rev. D 46, 510 (1992). 

[35] G. G. Raffelt and G. Sigl, "Self-induced decoherence in 
dense neutrino ga ses," Phys. Rev. D 75, 083002 (2007) 
[hep-ph/0701182| . 

[36] A. Esteban-Pretel, S. Pastor, R. Tomas, G. G. Raffelt 
and G. Sigl, "Decoherence in supernova neutrino trans- 
formations suppressed by deleptonization," Phys. Rev. D 
76, 125018 (2007) |arXiv:0706.2498l [astro-ph]]. 

[37] R. F. Sawyer, "The multi-angle instability in dense 
neutrino systems ," Phys. Rev. D 79 (2009) 105003 
|arXiv:0803.43T9l [astro-ph]] . 

[38] G. G. Raffelt, M. T. Keil, R. Buras, H. T. Janka and 
M. Rampp, "Supernova neutrinos: Flavor-dependent 
fluxes and spectra," lastro-ph/0303226 

[39] T. Fischer, S. C. Whitehouse, A. Mezzacappa et al, 



"Protoneutron star evolution and the neutrino driven 
wind in general relativistic neutrino radiation hydro- 
dynamics simulations," Astron. Astrophys. 517, A80 
(2010). |arXiv:0908.187T1 [astro-ph.HE]] . 

[40] L. Hudepohl, B. Muller, H. -T. Janka et al., "Neutrino 
Signal of Electron-Capture Supernovae from Core Col- 
lapse to Cooling," Phys. Rev. Lett. 104, 251101 (2010). 
|arXiv:0912.0"260l [astro-ph.SR]] . 

[41] J. F. Cherry, G. M. Fuller, J. Carlson, H. Duan and 
Y. Z. Qian, "Multi-Angle Simulation of Flavor Evolution 
in the Neutrino Neutronization Burst From an O-Ne- 
Mg Cor e-Collapse Super nova," Phys. Rev. D 82 (2010) 
085025 |arXiv:1006. 21751 [astro-ph.HE]] . 

[42] H. Duan and A. Friedland, "Self-induced suppression of 
collective neutrino oscillations in a supernova," Phys. 
Rev. Lett. 106, 091101 (2011) |arXiv: 1006.23591 [hep- 
ph]]. 

[43] G. G. Raffelt, "Self-induced parametric resonance in 
collective neutrino oscillations," Phys. Rev. D 78 (2008) 
125015 | [arXiv:0810.1407l [hep-ph]] . 

[44] Animated figures at http 

: //www. mppmu.mpg.de/supernova/ multisplits/ 

[45] B. Dasgupta and A. Dighe, "Collective three-flavor 
oscillations of supernova neutrinos," Phys. Rev. D 77, 
113002 (2008) |arXiv:0712.3798l [hep-ph]]. 

[46] M. T. Keil, G. G. Raffelt and H. T. Janka, "Monte 
Carlo study of supernova neutrino spectra formation," 
Astrophys. J. 590, 971 (2003) [astro-ph/0208035l . 

[47] B. Dasgupta, A. Dighe, A. Mirizzi and G. G. Raffelt, 
"Collective neutrino oscillations in non-spherical geome- 
try," Phys. Rev. D 78, 033014 (2008) |arXiv:0805.3300l 
[hep-ph]]. 

[48] H. Duan, G. M. Fuller, Y. -Z. Qian, "Symmetries in 
collective neutrino oscillations," J. Phys. G G36, 105003 
(2009). |arXiv:0808.2046l [astro-ph]]. 

[49] M. C. Gonzalez-Garcia, M. Maltoni and J. Salvado, 
"Updated global fit to three neutrino mixing: status 
of the hints of 8 1S > 0," JHEP 1004 (2010) 056 
[arXivTOOl. 45241 [hep-ph]]. 

[50] A. Esteban-Pretel, A. Mirizzi, S. Pastor, R. Tomas, 
G. G. Raffelt, P. D. Serpico and G. Sigl, "Role of dense 
matter in collective supernova neutrino transformations," 
Phys. Rev. D 78, 085012 (2008) [arXiv:0807.0659l [astro- 
ph]]. 

[51] B. Dasgupta, A. Dighe, A. Mirizzi and G. G. Raffelt, 
"Spectral split in prompt supernova neutrino burst: An- 
alytic t hree-flavor treatm ent," Phys. Rev. D77, 113007 
(2008). jarXiv:0801. 16601 [hep-ph]]. 

[52] G. Fogli, E. Lisi, A. Marrone and 1. Tamborra, 
"Supernova neutrino three-flavor evolution with dom- 
inant collective effects," JCAP 0904, 030 (2009). 
[arXiv:0812.3"03T1 [hep-ph]]. 

[53] H. Duan, A. Friedland, G. C. McLaughlin and R. Sur- 
man, "The influence of collective neutrino oscillations on 
a supernova r-process," arXiv: 1012.0532 [astro-ph.SR]. 

[54] |http://dollywood.itp.tuwien.ac.at/cernlib/| 



